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ABSTRACT 

A total of 235 Active Galactic Nuclei (AGN) from two different soft X-ray surveys 
(the ROSAT Deep Survey -DRS- and the ROSAT International X-ray Optical Survey 
-RIXOS-) with redshifts between and 3.5 are used to study the clustering of X- 
ray selected AGN and its evolution. A 2 sigma significant detection of clustering of 
such objects is found on scales < 40 — 80/i _1 Mpc in the RIXOS sample, while no 
clustering is detected on any scales in the DRS sample. Assuming a single power 
law model for the spatial correlation function (SCF), quantitative limits on the AGN 
clustering have been obtained: a comoving correlation length 1.5 < r$ < 3.3/i _1 Mpc 
is implied for comoving evolution, while 1.9 < r$ < 4.8 for stable clustering and 
2.2 < ro < 5.5 for linear evolution. These values are consistent with the correlation 
lengths and evolutions obtained for galaxy samples, but imply smaller amplitude or 
faster evolution than recent UV and optically selected AGN samples. We also constrain 
the ratio of bias parameters between X-ray selected AGN and IRAS galaxies to be 
< 1.7 on scales < 10 h~ x Mpc, a significantly smaller value than is inferred from local 
large-scale dynamical studies. 

Key words: surveys - galaxies: clusters: general - quasars: general - large-scale 
structure of Universe - X-rays: general 



1 INTRODUCTION 

Since the launch of ROSAT in 1991, a large number of sur- 
veys of soft X-ray selected sources have been undertaken 
using the Position Sensitive Proportional Counter (PSPC), 
with different sky coverages, depths and identification com- 
pletenesses. They range from the ROSAT All Sky Survey 
(RASS, Voges 1992) that covers the whole sky to a fairly 
bright flux limit, passing through the relatively large solid 
angle RIXOS (20 deg 2 Mason et al., in preparation), to the 
deeper ROSAT Deep Survey (henceforth DRS, Shanks et al. 
1991), Cambridge Cambridge ROSAT Serendipitous Survey 
(Boyle et al. 1995) or the UK Medium Survey (Carballo et 
al. 1995). The deepest and narrowest surveys, pushing to 
the limit the capabilities of the PSPC are the UK Deep Sur- 
vey (Branduardi-Raymont et al. 1994, McHardy et al. 1998) 
and the Deep Survey in the Lockman Hole (Hasinger et al. 
1993, Hasinger et al. 1998, Schmidt et al. 1998). 

Together, these surveys have given considerable insight 
into the evolution of the AGN X-ray luminosity function 



(Boyle et al. 1994, Page et al. 1996, Jones et al. 1997), their 
X-ray spectra (Bade et al. 1995, Carballo et al. 1995, Al- 
maini et al. 1996, Romero-Colmenero et al. 1996, Ciliegi et 
al. 1997, Mittaz et al. 1998) and their relation to Narrow 
Emission Line Galaxies (NELGs, Boyle et al. 1995, Griffiths 
et al. 1996, Page et al. 1997). 

However, so far the only direct study of the clustering 
properties of X-ray selected AGN is that by Boyle & Mo 
(1993, henceforth BM). They studied the local (z < 0.2) 
AGN in the Extended Einstein Medium Sensitivity Survey 
(EMSS, Stocke et al. 1991) and found no signal above 1 
sigma, consistent with the clustering properties of the UV 
and optically selected AGN. 

In this work we present for the first time a direct study 
of the clustering of soft X-ray selected AGN from ROSAT, 
using RIXOS and the ROSAT Deep Survey (positions and 
redshifts for the AGN identified in the five DRS fields have 
been kindly provided by O. Almaini, prior to publication). 

The samples are introduced in Section [| We discuss the 
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Figure 1. Histogram of the distribution of the redshifts of the 
AGN in RIXOS and DRS (see text). 

evidence for clustering in Section ^, using different model in- 
dependent methods. Quantitative measurements on the cor- 
relation strength and evolution are then obtained in Section 
^| These results are presented in Section and discussed in 
Section |^. We summarize our conclusions in Section 0. 

We have used Ho = 100 h km s^ 1 Mpc" 1 and go = 0.5 
(unless otherwise stated) throughout this paper. 



2 THE DATA 

We have used two large complete soft X-ray AGN samples: 
DRS (Boyle et al. 1994, Shanks et al., in preparation), and 
RIXOS (Page et al. 1996, Mason et al., in preparation). 

The DRS sample is the deepest, being a 'pencil beam' 
style survey in a few chosen directions in the sky, while 
RIXOS is wider (a collection of many shallower pencil 
beams), being a compromise between depth and surveyed 
solid angle. Both surveys together cover a wide range of red- 
shifts (see Fig. [ij) and allow an investigation of the clustering 
properties of X-ray AGN. 

DRS: Five deep PSPC exposures have been searched for 
sources down to the sensitivity limit of that instrument over 
the inner field of view (off-axis angle <18 arcmin). Due to 
the increasing width of the Point Spread Function (PSF) 
with off-axis radius and to vignetting, the flux limit of this 
survey is a function of off-axis radius and different for each 
field. This has been taken into account in our study, as well 
as the fact that the source identification completeness also 
varies with flux (see Section 3.2). A total of 107 X-ray AGN 
have been detected in these fields, with redshifts between 
and 3.5. 

RIXOS: This is a wider solid angle shallower survey than 
the DRS, in which the central 17 arcmin radius solid an- 
gle of 80 medium-deep PSPC pointings with | b \> 20° 
and exposure times longer than 8000 s have been source 
searched. About 90 per cent of the sources in a total of 
65 fields have been identified down to a uniform limit of 



Table 1. RIXOS fields used. The column labelled 'RIXOS' gives 
the RIXOS Field identification number and the column labelled 
'#' is the ROSAT sequence identification number (see Section 0). 
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Table 2. The different samples. 'TV' are the number of AGN 
in each sample. The column labelled '%' are the percentages of 
Poisson simulations with lik elihood values higher than those of 
each sample (see Section 3.1). The column labelled 'D c ' gives an 

distance between the objects in 



estimate of the mean comoving 

each sample (see Section |^) . The samples labelled with an asterisk 
are the ones used in Sections ^ and |i[ 

1 Mpc) 
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the sensitivity limit of all fields (Mason et al., in prepara- 
tion). We have further selected those fields with an exposure 
time longer than 10000 s, resulting in 43 fields (see Table [l]) 
and 128 AGN with redshifts between and 3.5. 

The median luminosities of the RIXOS and DRS sam- 
ples are similar Lo.5-2kcV ~ 0.2 x 10 44 ergs _1 (assuming a 
power law spectrum with an a — 1 energy index). The av- 
erage values are however different (L44) ~ 0.46 for RIXOS 
while (L44) ~ 0.23 for DRS. This indicates that, while the 
sources in both samples have similar overall luminosities, 
the luminous AGN in RIXOS are brighter than those in the 
DRS (see Figs. 2 of Boyle et al. 1994 and Page et al. 1996). 
We will se e th at this does not have any effect on our results 



in Section 3.2 



3x10 14 erg cm 2 s 1 (0.5-2 keV), significantly higher than 



The number of objects in each sample and the average 
redshifts for different redshift ranges are shown in Table H. 
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Also given are D c the comoving distances below which the 
total number of observed pairs (see Section fe.2| ) equals half 
the number of objects in the sample. This is an estimate of 
the mean distance between the objects in each sample, given 
the 'pencil beam' geometry of each of the ROSAT fields 
that make up our two surveys. D c is missing for the R1XOS 
1 < 2 < 2 sample because there are only 16 pairs in total in 
that sample. Our default samples are those marked with an 
asterisk in Table ^| the total RIXOS and DRS samples, a 
total of 235 AGN. These two samples have complementary 
redshift distributions: most RIXOS AGN are at z < 1 ({z) = 
0.84) while most DRS AGN are at z > 1 ((2) = 1.43) (see 
Fig. and Table |). In this sense, RIXOS reflects the 'local' 
behaviour of AGN clustering, while DRS can constrain its 
'medium-high redshift' evolution. 



bly due to a combination of the falling sensitivity (typical 
in a flux-limited survey) and the lower clustering ampli- 
tude at higher redshift ('positive evolution' see Section |J). 
This method does not detect any significant clustering in 
the whole RIXOS sample (0 < 2 < 3.5). 

This test has also been used to check for possible ef- 
fects of the galactic absorption and/or exposure times on 
the mean density of sources in different fields. We have found 
that there is not any significant difference between the low 
and high column density fields, nor between the shorter and 
longer exposure time fields, in terms of the surface density of 
RIXOS sources. This test was also repeated for the faint and 
bright sources separately, finding the same negative result. 
We can thus be confident in the uniformity of our source 
sampling with respect to 'instrumental' selection effects. 



3 IS THERE ANY EVIDENCE FOR 

CLUSTERING IN THESE SAMPLES? 

We have investigated the presence of clustering in the 
RIXOS sample in two different model-independent ways: 
by comparing the distribution of the number of sources per 
field in each sample with that expected from a purely Pois- 
son distribution (a variant of the counts-in-cells method, 
but in angular rather than spatial cells) and by comparing 
the total number of pairs of sources separated by a comov- 
ing distance r c with that expected from a uniform spatial 
distribution of sources. This second method has also been 
used for the DRS sample, the first not being adequate due 
to the small number of fields (5) of that sample. 



3.1 Counts— in— cells: checking the uniformity of 
our RIXOS sample 

If N is the total number of AGN in a sample, n is the to- 
tal number of fields (n = 43 for RIXOS), /j, = N/n is the 
observed average number of sources per field, and Ni is the 
number of sources in field i, we define a likelihood 



L = -^logP M (iV i ) 



(1) 



where P ll (Ni) is the Poisson probability of finding Ni from 
a Poisson process of average /j, 



Nl 



(2) 



For each sample, 1000 Poisson simulations with the 
same jj,, N and n as the real sample are performed, and 
the likelihood L is calculated for each one of them. The per- 
centage of Poisson simulations with likelihood values larger 
than those of the corresponding observed samples are given 
in the last column of Table |^. 

We can see that 96 per cent of the Poisson simulations 
have 'better' likelihood than the RIXOS 0.5 < 2 < 1 sample 
or, in other words, there is evidence for clustering at the 2 
sigma level in RIXOS in the 0.5-1 redshift interval. This is 
independent of the nature of the clustering. 

The significance becomes smaller if we consider together 
all AGN with < 2 < 1, probably because the volume sam- 
pled at low redshift is much smaller. The lack of cluster- 
ing signal in the higher redshift bin (1 < 2 < 2) is proba- 



3.2 Pairs of sources 

The previous method discards all the information on the 
spatial separation of the sources. It is obvious that two 
sources in the same field, but at the opposite ends of the 
redshift interval, are physically unrelated. The counts-in- 
cells method used above does not have a way of discriminat- 
ing against such cases, unless the redshift intervals are made 
smaller in which case the quality of the statistics worsens. 

We have performed new simulations in which as many 
sources as in the real samples are redistributed at random 
among the different fields (but not between the two different 
samples), keeping their redshifts and fluxes, but randomiz- 
ing their off-axis angles and 'azimuths' (angle between the 
line joining the field centre to the source and the meridian 
through the field centre). This has been done in a different 
way for the two samples: 

DRS: In this sample, the survey effective solid angle 
was different for every field and a function of off-axis an- 
gle within every field (see Table 2 of Boyle et al. 1994). 
A total of 25 different regions were defined, 5 for every 
field, corresponding to the off-axis ranges 0-10 arcmin, 10- 
12 arcmin, 12-14 arcmin, 14-16 arcmin and 16-18 arcmin. 
The flux limit Sum for source detection was different in 
each one of these 25 regions, as was the overall identifica- 
tion completeness for sources with S > 5n m . The 'effec- 
tive' solid angle f2 e ff(<Sii m ) °f the survey for S > Sii m is de- 
fined so that n c ff (Siim)/f2geom(Siim) i s the fraction of iden- 
tified sources down to that flux limit, where n goom (S'ii m ) 
is the total solid angle surveyed, also down to the same 
flux limit. When a simulated source with flux S' was ex- 
tracted, we found the highest Sy m to that S' > S[ im . The 
source could be detected over the total surveyed area with 
Shm < S[ im . A particular region within that area (and hence 
a field) was assigned at random to the source, proportion- 
ally to f2 ff(Sii m )/r2 c ff (S' liln ). Finally, the off-axis angle and 
azimuth were obtained assuming an uniform source distri- 
bution within the corresponding off-axis range. 

RIXOS: The source density is uniform within every field 
and field to field in this sample. Therefore, we assigned at 
random a RIXOS field to every source, and an off-axis angle 
(within 17 arcmin) and azimuth within that field. 

The number iV b s of pairs of sources in the two (sep- 
arate) real samples with comoving distance < r c is then 
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Figure 2. Complement to one of the significance of clustering as 
a function of the comoving distance r c , for the RIXOS (solid line) 
and the DRS AGN (dotted line). 1 and 2 sigma significance levels 
are shown with dashed lines 

obtained, as well as the mean number (among 100000 simu- 
lations) of random pairs Np up to the same separation. The 
distances r c are calculated in comoving coordinates accord- 
ing to the expressions given by Osmer (1981), for go = and 
0.5. Only pairs of sources in the same (real or simulated) field 
are used, to facilitate the calculation of the volume integrals 
introduced in Section ^j. 

We have used the integrated Poisson probability of ob- 
taining < A?ob s pairs for a distribution of mean y,{= Np): 

N abs -1 

P^(<N obs ) = J2 P "( f > ( 3 ) 

JV=0 

The complements to one of these probabilities (i.e. 
1 — Pp) are plotted in Fig. |^, for the r c values corresponding 
to each real pair in the RIXOS sample and the DRS sample. 
It is clear from this plot that there is a detection of clus- 
tering (at the ~95 per cent or 2 sigma level) in the RIXOS 
sample, at comoving distances r c < 40 — 80 ft -1 Mpc (e.g., 
iVobs = 22 pairs observed for Np = 14.6 pairs expected for 
r c < 83.03 h^ 1 Mpc). We have checked that the sources con- 
tributing to this signal are unrelated to the targets of the 
corresponding ROSAT pointings. As a general comment, the 
targets are mostly at z < 0.2 while the sources contributing 
to the signal have 0.2 < z < 1.4. 

It is also clear that there is not a significant detection 
of clustering in the DRS sample at any separation (e.g., for 
r c < 5.68 h' 1 Mpc, N ohs = 5 and Np = 4.85). None of these 
results changes significantly if go = is used instead of the 
default go = 0.5. 

We also note the lack of significant clustering at 
the smallest separation of any RIXOS pairs at (r c < 
6.44 ft" 1 Mpc, AU S = 1 and Np = 1.44). This fact, along 
with the non-detection in the DRS sample at similar scales 
and the detection of clustering at r c < 80 h~ Mpc will be 
used in the next section to constrain the clustering ampli- 
tude for different clustering evolution models. 

The lack of very close pairs is not due to the limited 



angular resolution of ROSAT. We have repeated the calcu- 
lation of both the real and simulated pairs excluding those 
pairs with an angular distance smaller than 1 arcmin (larger 
than the minimum distance at which the ROSAT-PSPC 
could resolve two separate sources at the flux levels relevant 
here), finding very similar results. 

We have investigated if the detection of clustering in 
RIXOS and not in the DRS could be due to the higher lu- 
minosity of the RIXOS sources. The average luminosity of 
the RIXOS sources contributing to the r c < 100/i _1 Mpc 
pairs is (L44) = 0.25, while that of the DRS sources at the 
same separations is (L44) = 0.24, the median luminosities 
being L44 = 0.17 and L44 = 0.19, respectively. Clearly, there 
is not any significant difference in luminosities between the 
two samples to which the observed difference in clustering 
strength could be attributed. However, there is a difference 
in their redshift distributions: all of the RIXOS sources at 
those separations are at z < 1.4, while most of the DRS 
sources have z > 1. This is not very surprising, considering 
that most models of structure formation predict the growth 
of inhomogeneities with cosmic time, or less clustering at 
higher redshift. 

Finally, we would like to emphasize that this test is 
sensitive to all pairs with comoving separations smaller or 
equal to the value of r c shown. It is also important to stress 
the independence of this test on any particular clustering 
model. 



4 THE INTEGRATED SPATIAL 
CORRELATION FUNCTION 

Since there is an excess of AGN pairs up to a certain co- 
moving separation, with respect to the expectation from a 
uniform distribution, we include now the possible spatial 
clustering of the sources. Indeed, the expected number of 
pairs up to a certain comoving separation will depend on 
the source correlation function integrated up to the corre- 
sponding comoving separation and taking into account the 
detailed geometry of the volumes sampled by the RIXOS 
and DRS observations. 

The spatial (3-dimensional) correlation function £(r) is 
defined as an excess probability (Peebles 1980). If n is the 
spatial density of the sources under study, the probability of 
finding a source in the volume dVi and another one in (IV2, 
separated by ri2, is: 

5P = n 2 (l + £(ri2))dVidy-2 (4) 

As mentioned above, the relevant quantity is the integrated 
spatial correlation function (SCF): 

f(0 = ± j v dV ^ r ) ( 5 ) 

where V is the (comoving) volume over which the pairs are 
counted. In the present paper, that volume is the intersec- 
tion of a sphere of radius r c with the cone defined by the 
maximum off-axis angles of each one of our ROSAT fields 
and the redshift limits of our samples. We call Vi to such 
volume around each of our sources with redshift Zi. 

For the DRS sample, this is further refined by slicing 
this sphere along the conical shells defined for the different 
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Figure 3. The 'best fit' ro (solid line) as a function of the co- 
moving separation r c for stable clustering for the RIXOS AGN, 
as well as the 2 sigma limits (dotted lines). 



effec tive solid angles and flux limits for every field (see Sec- 
tion 3.2), and summing up the integrals over those shells 
weighting them by 



, _ N{> 5 Um ) n e{t (Sii 

J — tt; - r X 



N(>S„ 



{Sun 



(6) 



s ), n eS (Sii m ) is 
geom(Sllm) is the to- 



where N(> S) is the surface density of sources with flux 
greater or equal to S, Su m is the flux limit of the corre- 
sponding off-axis shell, Smin is the lowest flux limit of the 
DRS (So.5-2keV = °-32 x fCT 14 ergcm 
the effective solid angle at Sum, and fi. 
tal surveyed solid angle at Sij m . This takes into account the 
fraction of AGN missed by the DRS due to the different flux 
limits at different off-axis radii, and the fraction missed due 
to the flux-dependent identification incompleteness of that 
survey. 

The expected mean number of pairs in a volume V is 
given by (Peebles f980): 



Ac. 



JVp(l + *(fc)) 



(7) 



Given a model for the SCF, the above integrals can be 
performed, and the parameters of the model constrained us- 
ing Pjv cxp (< A b s ) (see equation^). For example, a model 
with an amplitude too large would produce a value of A exp 
too high compared to A b s that would then be highly un- 
likely. 

We have assumed the usual power-law shape for the 

SCF 



£(r c ,z) = (l + z)- p (r c /r y 



(8) 



where ro is the correlation length, y = 1.8 (the results pre- 
sented below do not change significantly if 7 = 1.6 is used 
instead) and p is an evolutionary parameter (BM). If p = 
the clustering is constant in comoving coordinates ('comov- 
ing evolution'). If the clustering is constant in physical coor- 
dinates ('stable clustering') then p — 1.2. Linear growth of 
structures corresponds to p = 2 in an Q = 1 universe, higher 
values of p representing even faster non-linear growth. 



The evolution of the SCF is often parameterized, in 
physical coordinates, as 



C(r, 2 ) = (l + 2 )- (3+E) (r/r )- 



(9) 



where r is the physical distance and e = 7 + p — 3 (Groth & 
Peebles 1977). 
Hence: 



A cxp (p,r ) = ^ Np(Az)x 



x 1 + 



z;€Az 



(i + ^) (_p) 



Vi 



dVrZ 



(10) 



where Np(Az) is the number of poisson pairs within a red- 
shift interval Az (Np = X/Az Np(Az)) from the simula- 
tions. The sum ^ z is performed over all the sources with 
redshifts in Az, and it is essentially an average within Az 
of £(r c ) taking into account the clustering evolution model 
assumed in each case. The sum over Az is to deal with the 
changing spatial density of our sources as a function of red- 
shift, plus t he g eometry of our surveys and their flux limits 
(see Section 3.2). 

We show in Fig. ^| the 'best fit' ro (the value of ro that 
makes A cxp = A Q b s ) and the 2 sigma limits on ro for sta- 
ble clustering (p = 1.2), as a function of r c for the dis- 
tances between each real pair of RIXOS AGN (each new 
point in the lines uses the cumulative distribution of pairs 
up to distance r c ). We can see that there are only upper 
limits to the value of ro at small separations, but the re- 
quired value of ro becomes different from zero at separa- 
tions r c < 40 — 80 h~ Mpc, in agreement with the findings 



of Section 3.2 



Similarly, each pair (p,ro) can be assigned a probability 
for a fixed value of r c using equations |^ and ^| This will be 
used in the next section to constrain the clustering of X-ray 
selected AGN and/or its evolution. 

In Figs. 2 and 3 we see that we have to go up to co- 
moving separations of ~ 40 ft -1 Mpc to find a significant 
excess of pairs with respect to an unclustered population of 
sources. Fig. 3 also shows that the required comoving corre- 
lation length is much smaller, so the detection of the cluster- 
ing signal occurs mostly at the tails of the correlation func- 
tion where clustering is weak. We do not detect an excess 
of pairs at comoving separations r c < 10 h^ 1 Mpc probably 
because of the relatively low density of objects in our sam- 
ple, in which case the number of expected pairs is always 
small at small separations. Although all this might seem 
puzzling, Fig. 3 shows that a comoving correlation length 
ro ~ 5 h^ 1 Mpc is consistent at all separations. 



5 LIMITS ON THE CLUSTERING STRENGTH 
FOR DIFFERENT EVOLUTION MODELS 

We show in Fig. |j| the 1, 2 and 3 sigma contours in the 
(p,ro) space for the distance of the closest observed pair in 
the RIXOS sample (r c < 6.44 h -1 Mpc). Since no clustering 
is detected at those separations, only upper limits to the 
value of ro for each p can be obtained, ranging from ro < 
3.8 Mpc for comoving clustering to ro < 5.5 Mpc for 
linear growth (2 sigma). 



© 0000 RAS, MNRAS 000, 000-000 



6 F.J. Carrera, X. Barcons, A.C. Fabian et al. 



Figure 4. Probability contour levels corresponding to 1 sigma 
(solid line), 2 sigma (dashed line) and 3 sigma (dot-dashed line) 
in the (p, ro) space for the RIXOS AGN at distances r c < 
6.44 h' 1 Mpc. 



Figure 6. Probability contour levels corresponding to 1 sigma 
(solid line), 2 sigma (dashed line) and 3 sigma (dashed-dotted 
line) in the (p, ro) space for the RIXOS AGN at distances r c < 
83.03 h' 1 Mpc. 



-4 




Figure 5. Probability contour levels corresponding to 1 sigma 
(solid line), 2 sigma (dashed line) and 3 sigma (dashed-dotted 
line) in the (p,ro) space for the DRS AGN at distances r c < 
5.68 h' 1 Mpc. 



The 1, 2 and 3 sigma contours are plotted in Fig. ^ 
for the DRS AGN at r c < 5.68 h" 1 Mpc. This distance has 
been chosen to represent similar scales to those in the first 
RIXOS pair, but, due to the absence of any signal of clus- 
tering in the DRS sample, it is representative of the whole 
set of distances. Only upper limits are found. They are more 
restrictive than those from RIXOS for p < 0, because DRS 
is deeper than RIXOS. Hence, it can constraint better the 
high redshift behaviour of clustering (p < implies stronger 
clustering at higher redshift). 

The allowed region of the (p, ro) space from the clus- 
tering signal in RIXOS is shown in Fig. |^. The comoving 
distance at which this is calculated (r c < 83.03 /t _1 Mpc) is 



Figure 7. Probability contour levels corresponding to 2 sigma 
from RIXOS r c < 6. 44 h' 1 Mpc (solid line) and r c < 
83.03/i _1 Mpc (dashed lines), and DRS r c < 5.68h~ 1 Mpc 
(dashed— dotted line). The allowed region is that below the solid 
and dashed— dotted lines and between the two dashed lines (shown 
shaded). 



the largest one below which there is any clustering signal (at 
2 sigma) in RIXOS, and is taken as a representative value, 
ro = (no clustering) is excluded at >2 sigma level for any 
value of p, in accordance with our detection of clustering. 

Figs. ^, |B| and ^ can be combined to produce a joint set 
of limits on the clustering strength and evolution, using the 
RIXOS small and large scale 'low redshift' results and the 
DRS small scale 'high redshift' results. The 2 sigma contours 
from those three plots are shown in Fig. [?] together. Note 
that a region in the (p, ro) parameter space 'outside' any 
of the 2 sigma contours is actually excluded at more than 
2 sigma significance, when the three sets of contours are 
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taken into acount simultaneously (or any combination of 
two of them). 

To reproduce simultaneously the lack of close pairs and 
the clustering below ~ 80 ft -1 Mpc in RIXOS, ro has to be 
between the lower dashed line and the solid line in Fig. ^ 
This translates to 1.5 < ro < 3.8 /i" 1 Mpc for comoving 
evolution, 1.9 < ro < 4.8 for stable clustering and 2.2 < ro < 
5.5 for linear evolution of clustering (>2 sigma significance). 

The addition of the DRS limits reduces slightly the up- 
per limit on ro for the comoving evolution model: 1.5 < ro < 
3.3 Mpc, leaving the others unchanged. 

Using go = increases the limits: comoving evolution 
is allowed for 2 < r < 4.2 h^ 1 Mpc, 2.5 < r < 5.8 for 



stable clustering and 2.9 



< 



ro 



< 8.3 for linear evolution of 



clustering. These are also >2 sigma limits. 

We have addressed the question of how much our results 
depend on the exact r c value chosen to reflect the scales at 
which there is some signal in our RIXOS sample. We have re- 
peated the calculations for r c < 43.96 /i" 1 Mpc, the first co- 
moving distance at which there is any signal above 2 sigma. 
They are very similar to those for r c < 83.03 ft -1 Mpc, the 
lower 2 sigma limits being ~19 per cent higher. We there- 
fore conclude that r c < 83.03 ft" 1 Mpc is representative of 
the whole r c < 40 - 80 h' 1 Mpc range. 



6 DISCUSSION 

The clustering of galaxies seems to evolve according to the 
stable model or even faster with correlation lengths ro ~ 
3 — 7 Mpc. Smaller correlation lengths (ro < 3 ft -1 Mpc) 
would be required if the evolution is comoving (see e.g. In- 
fante, de Mello & Menanteau 1996, Le Fevre et al. 1996, 
Hudon & Lilly 1996, Carlberg et al. 1997, Brainerd & Smail, 
1998). These values are of the order of our limits for the cor- 
responding evolution models, and therefore our results imply 
similar clustering properties of galaxies and X-ray selected 
AGN. 

However, clustering of optically-UV-selected AGN ap- 
pears to be stronger than that implied by our results. For in- 
stance, Croom & Shanks (1996) find r = 5.4 ± 1.1 h' 1 Mpc 
and comoving evolution assuming a 7 = 1.8 power law. Us- 
ing a biasing model they obtain ro ~ 7 — 8/i _1 Mpc and 
slow evolution (comoving or slightly faster). These results 
were obtained at r c < 10 ft -1 Mpc. They are higher than 
our RIXOS and DRS upper limits at comparable scales. 

A number of recent works have found an increase in 
the clustering of AGN for increasing redshift: Stephens et 
al. (1997) found ro = 18 ± 8/i _1 Mpc assuming comoving 
clustering evolution with a sample of 56 z > 2.7 AGN over 
~22 deg 2 . This is much higher than our corresponding up- 
per limits, and assuming evolution in comoving coordinates 
would lead to p < 0, again in conflict with our results. How- 
ever, we only have 5 AGN with z > 2.7, so we essentially 
don't have any information at these redshifts. A cluster- 
ing evolution scenario in which clustering is strong at high 
redshifts, then it decays, and grows again at lower redshift 
would in principle be compatible both with their and our 
results. This is qualitatively the behaviour of the model 
suggested by Bagla (1998), in which higher (rarer) mass 
overdensities collapse early and cluster very strongly. The 
clustering amplitude then decreases while lower and lower 



mass objects collapse. When the average mass objects have 
collapsed, clustering starts growing again because of their 
mutual gravitational attraction. 

La Franca, Andreani & Cristiani (1998) found a 2 sigma 
significant increase in the quasar (AGN with A4b < 23) clus- 
tering amplitude between z — 0.95 and z = 1.8 using a new 
sample of objects. However, if other samples are also taken 
into account, the significance of this decreases. In partic- 
ular, if we use the four values for £(r c < 15 h~ Mpc, z) 
that are mutually independent in their work and fit £(15, z) 
to those data, p > cannot be excluded at more than 75 
per cent probability. The data points used in this fit are: 
f(15,0.97) = 0.5 ± 0.2 and f(15, 1.85) = 0.8 ± 0.3 from 
La Franca et al. (1998), £(15,0.05) = 0.2 ± 0.3 from Boyle 
& Mo (1993) and Georgantopoulos & Shanks (1994), and 
£(15,3.1) = 1.2 ± 0.7 from Kundic (1997). Moreover, a 
2 sigma effect as the one reported by La Franca et al. (and 
indeed, as our own clustering detection) is not very signifi- 
cant and requires further work to confirm or reject it. 

An interesting conclusion of our results applies to the 
estimates of the anisotropics introduced in the X-ray back- 
ground (XRB) by source clustering. Several studies have 
used the angular Auto Correlation-Function (ACF) of the 
X-ray Background to constrain the contribution of AGN 
to the XRB (Carrera & Barcons 1992, Carrera et al. 1993, 
Georgantopoulos et al. 1993, Danese et al. 1993, Chen et 
al. 1994, Soltan & Hasinger 1994). In general, those stud- 
ies coincided in stating that sources clustered on scales of 
6 to 8/i _1 Mpc with a correlation fixed in comoving coordi- 
nates could not produce more than about half of the XRB. 
A population of sources with a smaller correlation length or 
faster evolution (stable or linear) could make up the remain- 
ing XRB. These constraints are relaxed by our work, since 
it appears that X-ray selected AGN present weaker cluster- 
ing and a faster than comoving evolution in their correlation 
function. 

Recent soft X-ray surveys show that soft broad-line 
AGN only contribute ~50-60 per cent of the soft XRB, and 
that the contribution from harder sources (NELGs or ab- 
sorbed AGN) grows at faint fluxes (Boyle et al. 1994, Page 
et al. 1996, Romero-Colmenero et al. 1996, Almaini et al. 
1996). Whatever the ultimate nature of the hard sources 
turns out to be, if their clustering properties are like those of 
galaxies or like those of the AGN studied here, the whole soft 
XRB intensity could be produced by AGN and NELGs with 
stable or linear clustering evolution and ro values within the 
limits found here (see e.g. Fig. 5 of Carrera & Barcons 1992). 

It is also interesting to assess the impact of our studies 
on the contribution of AGN to the hard XRB. If hard X-ray 
AGN cluster as the soft X-ray AGN considered here, they 
could produce the whole hard XRB without exceeding the 
upper limits from the ACF (Carrera et al. 1993, Danese et 
al. 1993). At the moment, hard X-ray surveys only resolve 
a small fraction of the hard XRB (Ueda et al. 1998) . 

Finally, we can use our upper limits on ro from the 
smallest separation RIXOS pair to limit the ratio of bias 
parameters between X-ray selected AGN and IRAS galaxies: 
bx/bi ~ cr&x/a&i, where: 



2 

C-8 



72(r /8r 



(11) 



27(3- 7 )(4- 7 )(6-7) 
is the variance of the counts in a sphere of radius 8ft" 1 Mpc 
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for a power-law clustering model, and <jsi = 0.69 ± 0.04 
(Fisher et al. 1994). Our results for the three evolution 
models considered span bx/bi < 0.8 — 1.7. These up- 
per limits are somewhat smaller than previous estimates: 
bx ~ (6.8 ± l.6)Q ' 6 from the comparison of the dipole 
of bright 2-10 keV selected AGN with the motion of the 
Local Group (Miyaji 1994), bx/bi ~ 1.5 — 5 from cross- 
correlation of those AGN with IRAS galaxies (Miyaji 1994), 
and bx ~ 5.6 from a comparison between the X-ray and mi- 
crowave backgrounds (Boughn, Crittenden & Turok 1998). 
This difference might be due to the different scales sam- 
pled. Our result applies to scales ~ 10 /i -1 Mpc which is 
where fluctuations are usually normalised, whilst the other 
studies measure the bias parameter on much larger scales 
~ 1000 h~ Mpc. Our results are however similar to those 
of Treyer et al. (1998) {bx ~ 0.9 - 1.8) obtained from a 
study of the harmonic coefficients of the large angular scale 
fluctuations of the XRB from HEA01 A2 data. 



7 CONCLUSIONS 

The spatial correlation function of soft X-ray selected AGN 
has been studied here using data from two different X-ray 
surveys: the ROSAT Deep Survey (DRS) and RIXOS. 

Some indication of clustering at the 2 sigma level in 
RIXOS in the z = 0.5—1 redshift range has been found using 
a variant of the counts-in-cells method. A more powerful 
test has also been performed, by comparing the number of 
pairs of sources in each sample with the number of pairs 
expected from a uniform distribution of sources, finding a 
~2 sigma significant clustering signal in the whole RIXOS 
sample ((2) = 0.53) in comoving scales r c < 40—80 h" 1 Mpc. 
No significant detection has been found in the DRS sample 
((z) = 1.4) or in the RIXOS sample at small scales. Both 
tests are model-independent. 

Quantitative measurements of the clustering and its 
evolution have been obtained from the integrated spatial 
correlation function (essentially the number of pairs with 
comoving separations smaller or equal than r c ), assuming 
Poisson statistics, and a power-law shape for the correla- 
tion function, with slope 7 = 1.8 and correlation length ro. 

Combining the limits from RIXOS (r c < 6 Mpc and 
r c < 80 hT 1 Mpc) with the DRS limits (r c < 6 h' 1 Mpc), we 
obtain: 1.5 < r < 3.3 ft -1 Mpc for comoving clustering, 
1-9 < ro < 4.8 for stable clustering and 2.2 < ro < 5.5 
for linear evolution of clustering. Using go = we obtain 
~20-30 per cent higher limits. 

These results are compatible with the clustering prop- 
erties of 'normal' galaxies, but would imply weaker cluster- 
ing in X-ray selected AGN than in optically-UV-selected 
ones. Our results do not support an increase in the cluster- 
ing amplitude with redshift. However, we do not have many 
sources above z > 2.5, so we cannot rule out neither support 
the strong clustering above that redshift, found by Stephens 
et al. (1997). 

AGN (and/or NELGs) clustered like the sources studied 
here could produce most of the hard and soft X-ray back- 
ground without exceeding the observed limits on its Auto 
Correlation Function. 

The lack of very close pairs in RIXOS implies that the 
ratio of bias parameters between X-ray selected AGN and 



IRAS galaxies (where there is evidence that biasing is small, 
bi w 1) is bx/bi < 0.8 — 1.7, somewhat smaller than previous 
results. 

A 3 sigma detection at comoving separations r c < 
lO/i^Mpc would necessitate ~ 100 - 110 'RlXOS-like' 
fields for comoving evolution and ro = 6/i _1 Mpc. This is 
especially relevant for future X-ray survey identification pro- 
grammmes (such as the XMM-XID) with a view to detect 
AGN clustering. It is much more efficient to use a com- 
pact connected area (circular or elongated) for this purpose, 
rather than serendipitous pointings distributed all over the 
sky. In this sense, the ROSAT All Sky Survey identification 
programmes will be very useful to provide a local 'anchor' 
for X-ray AGN clustering studies (e.g. Zickgraf et al., 1998). 
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